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Abstract 

We show that, in the low-scale type-I seesaw model, renormalization group running of neutrino 
parameters may lead to significant modifications of the leptonic mixing angles in view of so-called 
seesaw threshold effects. Especially, we derive analytical formulas for radiative corrections to 
neutrino parameters in crossing the different seesaw thresholds, and show that there may exist 
enhancement factors efficiently boosting the renormalization group running of the leptonic mixing 
angles. We find that, as a result of the seesaw threshold corrections to the leptonic mixing angles, 
various flavor symmetric mixing patterns (e.g., bi-maximal and tri-bimaximal mixing patterns) can 
be easily accommodated at relatively low energy scales, which is well within the reach of running 
and forthcoming experiments (e.g., the LHC). 
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I. INTRODUCTION 



Experiments on neutrino oscillations have opened up a new window in searching for 
new physics beyond the Standard Model (SM) during the past decade. Since neutrinos are 
massless particles within the renormalizable SM, one often extends the SM particle content 
in order to accommodate massive neutrinos. Among various theories giving rise to neutrino 
masses, the seesaw mechanism attracts a lot of attention in virtue of its natural and simple 

nil 

description of tiny neutrino masses. In the conventional type-I seesaw model [lH4(, three 
right-handed neutrinos are introduced with super-heavy Majorana masses far away from the 
electroweak scale A EW = (9(100) GeV. Small neutrino masses are strongly suppressed by 
the ratio between the electroweak scale and the large mass of the right-handed neutrinos, 
which on the other hand, leaves the theory lacking in experimental testability. However, 
there are alternatives that allow us to realize the seesaw mechanism at an experimentally 
accessible level, e.g., the TeV scale. One popular way to lower the seesaw scale is to introduce 
additional singlet fermions, which have the same masses as the right-handed neutrinos but 
different CP parity They can be combined with the right-handed neutrinos to form 
four- component Dirac fields, while the lepton number is broken by a small Majorana mass 
insertion. In such a scenario, the masses of the light neutrinos are strongly suppressed by the 
small Majorana mass insertion instead of the seesaw scale, while the non-unitarity effects in 
the lepton flavor mixing can also be boosted to an observable level Other possibilities 



like structural cancellation 







lOj or the minimal flavor seesaw mechanism 
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121 ] may be 



employed to construct low-scale seesaw models. 

In principle, the current neutrino parameters are observed via low-energy neutrino os- 
cillation experiments. On the other hand, the seesaw- induced neutrino mass operator is 
usually given at the seesaw scale. Therefore, neutrino parameters are subject to radiative 
corrections, i.e., they are modified by renormalization group (RG) running effects. Typically, 
at energy scales lower than the seesaw threshold (i.e., the mass scale of the corresponding 
seesaw particle), the RG running behavior of neutrino masses and leptonic mixing param- 
eters should be described in an effective theory, which is essentially the same for different 
seesaw models. However, at energy scales higher than the seesaw threshold, full seesaw 
theories have to be considered, and the interplay between the heavy and light sectors could 
make the RG running effects particularly different compared to those in the effective the- 
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ory. However, in the spirit of some grand unified theories (GUTs), a unified description of 
fermion masses and flavor mixing depends on the lepton flavor structure at the GUT scale, 
which inevitably requires the RG running between seesaw particle thresholds and above. 
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The full sets of renormalization group equations (RGEs) in the type- 1 lri4], type-II 
and type-Ill 19|] seesaw models have been derived, both in the SM, and in the Minimal 
Supersymmetric Standard Model (MSSM) 2CH26I]. The general feature of the running pa- 
rameters have also been intensively studied in the literature, and it has been shown that 
there could be sizable radiative corrections to the leptonic mixing parameters at superhigh 



energy scales (see, e.g., Refs. 27], |28| and references therein). In particular, certain flavor 



symmetric mixing patterns can be achieved at the GUT scale indicating that there might ex- 
ist some flavor symmetries similar to the gauge symmetry (see, e.g., Ref. 29( and references 
therein) . 

Since the RG evolution together with the threshold effects of the heavy seesaw particles 
may result in visible corrections to the leptonic mixing parameters, there is a need to look 
into the RG running effects on the leptonic mixing parameters in the TeV seesaw model. 
In this work, we will explore in detail the RG evolution of neutrino masses and leptonic 
mixing parameters in the TeV type-I seesaw model. In particular, we will show that the 
threshold effects play a crucial role in the RG running of neutrino parameters, and some 
phenomenologically interesting flavor symmetric mixing patterns may be feasible even at an 
observable energy scale. 

The remaining part of this work is organized as follows. In Sec. HIl we first present general 
RGEs of the neutrino parameters in the type-I seesaw model; in particular, the matching 
conditions in crossing the seesaw thresholds. Then, the threshold corrections to the neutrino 
parameters are discussed in detail in Sec. IIHI Based on the analytical results, we illustrate 
a phenomenologically interesting numerical example in Sec. IIV} in which the bi-maximal 
leptonic mixing pattern at TeV scale is shown to be compatible with experimental data 
when the seesaw threshold corrections are properly taken into account. Finally, a brief 
summary and our conclusions are given in Sec. |Vj In addition, there are two appendices 
including the detailed analytical treatment of the threshold effects. 
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II. RGES IN THE TYPE-I SEESAW MODEL 



The simplest type-I seesaw model is constructed by extending the SM particle content 
with three right-handed neutrinos u R = (vri, z/r 2 , ^3) together with a Majorana mass term 
of right-handed neutrinos. The mass part of the lepton sector Lagrangian reads 

-£ = h<pY e e R + 40}>r + 2^ M R^R + h - c -> W 

where <fi is the SM Higgs boson, while £^ and denote lepton doublets and right-handed 
charged leptons, respectively. Here Y e and Y v are the corresponding Yukawa coupling matri- 
ces with Mr being the Majorana mass matrix of the right-handed neutrinos. Without loss of 
generality, one can always perform a basis transformation and work in the basis in which Mr 
is diagonal, i.e., Mr = diag(M 1; M 2 , M 3 ). At energy scales below the right-handed neutrino 
masses, the heavy right-handed neutrinos should be integrated out from the theory and the 
masses of the light neutrinos are effectively given by a non-renormalizable dimension-five 
operator 

-£f 5 = ^(40)-M0 T ^)+h.c., (2) 
where the effective coupling matrix k can be obtained from the type-I seesaw formula 

k = Y u M^Yj, (3) 
and the Majorana mass matrix of the left-handed neutrinos is 

m u = kv 2 , (4) 

with v ~ 174 GeV being the vacuum expectation value of the Higgs field. 

Because of the seesaw thresholds, the RG running of neutrino parameters needs to be 
treated separately. At energies above the seesaw thresholds, the full theory should be con- 



sidered and the relevant beta-functions are given by 2CH23| 



(\Y 

167r 2 /i^ = (a e + C e e H e + C:H u )Y ej (5) 



.2, an 



^ = K + C e u H e + C V V H V ) Y v , (6) 



167r 2 /i^ = C R M R (YX) + Oi {YjY u ) T M R , (7) 



where H f = Y f Yj for / = e,v,u,d and the coefficients (C e e , C v e , C% C v v , C R ) 
(3/2, —3/2, —3/2, 3/2, 1) in the SM. The coefficient a u is flavor blind and reads 

a v = tr (3H U + 3H d + H e + H v ) - ( ^-gf + \g\ J . 



If we make use of m v at energy scales above the seesaw threshold, we can derive from Eqs. ([2D 
and (ED- © 

= m„ = 2a„m„ + (C*if e + C m H v ) m v + {C e v H e + C m ^) T , (9) 

with C m = 1/2. Here, for simplicity, we have defined t = ln/x/(167r 2 ). 

At energies below the seesaw thresholds, i.e., in the effective theory, neutrino masses are 
attached to the dimension-five operator, and the RGE of k reads 

k = a K K + (C e u H e ) k + k (C e u H e ) T , (10) 

where 

a K = 2tr (3H U + 3H d + H e ) + A - 3g 2 2 , (11) 

with A denoting the SM Higgs self-coupling constant. 

In crossing the seesaw thresholds, one should ensure that the full and effective theories 
give identical predictions for physical quantities at low energy scales, and therefore, the 
physical parameters of both theories have to be related to each other. In the case of the 
neutrino mass matrix, this means relations between the effective coupling matrix k and the 
parameters Y v and Mr of the full theory. This is technically called matching between the 
full and effective theories. For the simplest case, if the mass spectrum of the heavy singlets 
is degenerate, namely M\ = M 2 = M 3 = M , one can simply make use of the tree- level 
matching condition at the scale \x = M 

k\ Mo = Y v M R Yj\ Mo . (12) 

In the most general case with non-degenerate heavy singlets, i.e., M\ < M 2 < M3, the 
situation becomes more complicated and the heavy singlets have to be sequentially decou- 



pled from the theory [30|. For instance, at energy scales between the n-th and (n — l)-th 
thresholds, the heavy singlets are partially integrated out, leaving only a 3 x (n — 1) sub- 
matrix in Y u , which is non-vanishing in the basis, where the heavy singlet mass matrix is 




FIG. 1: Feynman diagrams of the one-loop A corrections to the neutrino mass operator in the 
effective theory (left) and in the full theory (right), where p, q and k are the corresponding momenta. 



diagonal. The decoupling of the n-th heavy singlet leads to the appearance of an effective 
dimension- five operator similar to that in Eq. fl2]), and the effective neutrino mass matrix 
below M n is described by two parts 



w!f = v 2 



(13) 



where (n) labels the quantities relevant for the effective theory between the n-th and (n — 1)- 
th thresholds. In the SM, the RGEs for the two terms above have different coefficients for 
the gauge coupling and Higgs self-coupling contributions, which can be traced back to the 
decoupling of the right-handed neutrinos from the full theory. For instance, we show in 
Fig. [1] a comparison between the A corrections in the effective theory (left diagram) and in 
the full theory (right diagram). The loop integral of the left diagram, i.e., 

d 4 k 1 1 



(27r) 4 k 2 — m 2 ^ (k — p — q) 2 — m 2 ^ 
is UV divergent, whereas the loop integral of the right diagram, i.e., 

d 4 k k + M n 1 1 



(2tt) 4 k 2 -M 2 (k + p) 2 - ml (q - k) 2 - m 2 
is UV finite. Therefore, there is no A correction to the neutrino mass operator in the full 



theory, while A enters the beta-function of k in the effective theory. See also Ref. [3l| for a 
detailed discussion. We remark that, if the relevant seesaw threshold is above the soft SUSY- 
breaking scale, such a mismatch is absent in the MSSM due to the supersymmetric structure 
of the MSSM Higgs and gauge sectors. Therefore, this feature may result in significant RG 
running effects only in the SM, in particular when the mass spectrum of the heavy neutrinos 
is quite hierarchical. 
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III. THRESHOLD EFFECTS IN RGES 



We continue to determine the threshold corrections to the neutrino parameters. Accord- 
ing to Eq. (TO?]) the neutrino mass matrix between the seesaw thresholds consists of two parts 
k and Y v M^Yj . The beta-functions of these two parts have different coefficients in the 



terms proportional to the gauge couplings and A 3l|- In the framework of the SM, flavor 
non-trivial parts in the beta- functions of k and Y v are dominated by H„ and H e , which are 
essentially negligible in the RG running. For example, in a basis where H e is diagonal, the 
largest entry of H e is proportional to the square of the tau Yukawa coupling, i.e., y^, which is 
at the order of 10~ 4 and too small to be significant in the RG running. As for H u , stringent 
constraints from the unitarity of the leptonic mixing matrix indicate that their contributions 
in the RG running are not comparable to the flavor blind parts, i.e., the gauge couplings 
and A. Hence, we neglect the flavor non-trivial parts in the beta-functions in the following 
analytical analysis. 

For the RG running of the neutrino parameters from one seesaw threshold M n to the 
nearby one M n _ 1; the two parts of the neutrino mass matrix are basically re-scaled to an 
and bY u M^Yj, with a ^ b. Explicitly, the mass matrix of the light neutrinos at fi = M n _i 
approximates to 



m,. 



± ~ bv 2 (k + Y U M R X T ) \ Mn + (a- b)v\\ Mn 



Til 



v \M n 



(afe" 1 - 1) 



2 

V K 



blm v \ Mn + ev k\ m ,, 



(14) 



where 



exp 



exp 



exp 



lnM„ 



1 hM„ 

1671"^ 




lnM„ 



lnM„ 



InM„ 



111 A I n 



a K dt 



2a u dt 



M 



n-1 



2a v 

KiTT 2 



[a. 



2a u ) dt 



M, 



n-1 



M„ 



16tt 2 



- 1. 



(15) 
(16) 
(17) 



In the limit M n _i = M n , one can easily observe that the second term of Eq. (|14j) disappears. 

In the low-scale seesaw model, when the evolution of the neutrino mass matrix crosses 
two thresholds (e.g., from M n to M n _x), the threshold corrections, i.e., the second term in 
Eq. (Jill), provides corrections to the neutrino mass matrix. In comparison to the RG running 



corrections, the threshold effects are much more significant in changing the value of the 
neutrino parameters. In order to see this point, we recall that when the neutrino parameters 
are evaluated below the seesaw scale, i.e., in the effective theory, the only flavor non-trivial 
corrections in the beta- functions are related to the charged lepton Yukawa coupling Y e . 
As we mentioned before, the largest entry in Y e is constrained by the tau mass, and can 
only be at percentage level. Therefore, no visible RG running effects can be gained in the 
effective theory. At energies above the seesaw scale, Y v contributes to the RG running of 
the neutrino mass matrix as shown in Eq. (Q. In the ordinary seesaw models with GUT 
scale right-handed neutrino masses, Y v can be sizable, e.g., at order unity. However, in the 
low-scale seesaw model, Y v should be relatively small in order to keep the stability of the 
masses of the light neutrinos. Although there exist mechanisms that could stabilize neutrino 
masses without the requirement of a tiny Y v , there is still very strong restrictions from the 
unitarity of the leptonic mixing matrix. For example, if one of the right-handed neutrinos is 
located around the electroweak scale Aew> the unitarity of the leptonic mixing matrix sets 
a general upper bound H v < 1CT 3 unless severe fine-tuning is involved. Therefore, large RG 
running effects on the mixing parameters can hardly be obtained above the seesaw scale. In 
the case of threshold corrections, the situation becomes quite different, since in Eq. (I14p the 
flavor non-trivial corrections are proportional to the gauge couplings and A, which do not 
suffer from the unitarity constraints and can be of order unity. In particular, if the light 
neutrino mass spectrum is highly degenerate, sizable corrections to the neutrino parameters 
could be naturally obtained. 

In order to investigate the threshold effects on the neutrino parameters, we present series 
expansions of the leptonic mixing angles and the light neutrino masses. In particular, we 



make use of the so-called Casas-Ibarra parametrization 



32j, in which Y v is determined by 



y v = u<Jd~ k o<Jm^, (18) 

in a basis where Mr is diagonal. Here, U is the leptonic mixing matrix D K denotes the 
eigenvalue matrix of k, while O = ^23 (^1)^13(^2)^12(^3) with Rij({} k ) being the elemen- 
tary rotations in the 23, 13, and 12 planes, respectively. Different from the leptonic mixing 
angles, $j are in general complex. 1 For simplicity, we will assume CP conservation in this in- 



Sec Rcf. 



in the MSSM. 



33j for a discussion on the role of the matrix O in the RG evolution of the neutrino parameters 
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vestigation (except from in the appendices), i.e., O is a real orthogonal matrix parametrized 
by three real angles. A detailed derivation with CP-violating effects can be found in Ap- 
pendix |A] The neutrino mass matrix at the scale /i = M n can be diagonalized by using a 
unitary transformation, viz. 

m v\ Mn = U diag(mi, m 2 ,m 3 ) U T , (19) 

with rrii {i = 1,2,3) being the light neutrino masses. Then, at the scale \x = M n _i, the 
first-order correction to the leptonic mixing matrix is calculated by 

{U T KU) kj _ ^ y/mym 



U aj = U aj + ev 2 V yu ^- U ak = U aj + e V ^J^O kn O jn U ak (20) 

and the masses of the light neutrinos are given by 

rhi = bnii + ev 2 (U T KU)u = (b + e0 2 in ) m*. (21) 

Here n denotes the radiative corrections between the n-th and (n — l)-th thresholds. See 
Appendix [A] for the detailed derivation of Eqs. (120]) and ff2T]) . 

Inserting the parametrizations of U and O into the above equations, one can easily obtain 
the threshold corrections to the leptonic mixing parameters. For example, evolving the RGEs 
from M3 to M 2 , the leptonic mixing angles are modified to 



■- Ois + ec lC2 ^T 3 1 2 12V + 12V , 22 
m 3 — m 2 m 3 — m\ 



6 12 ~ 12 + esis 2 c 2 - , (23) 

m 2 — mi 

/ x/tui x/rfki \ 

V23 6*23 + eci^y/rrvi -s 2 si 2 h SiC 2 Ci 2 , (24) 

\ m 3 — mi m.3 — m 2 / 

where Sj = sin-$j and q = cos^j, and the terms proportional to esi 3 have been ignored. The 
light neutrino masses at the scale \i = M 2 are approximately given by 

fh\ = (fo + £s 2 )mi, (25) 
fh 2 = (b + es^cl) m 2 , (26) 
m 3 = (6 + ec\c 2 2 ) m 3 . (27) 

We refer the reader to Appendix [B] for the complete analytical derivation of the threshold 
corrections to the leptonic mixing angles and the light neutrino masses. 
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Note that the mass differences in the denominator of Eqs. (j22l) - (!24l) could strongly enhance 
the threshold corrections. We give a rough estimate of the enhancement effects. Unless 
in some fine-tuning limits, kv 2 should generally be the same order of magnitude as the 
neutrino masses, i.e., kv 2 ~ m u . Taking into account the neutrino mass-squared differences 
Am 2 j = m 2 — m 2 measured by neutrino oscillation experiments, one can calculate that, 
if mi ~ 0.05 eV, kv 2 j{m\ — 7712) = (9(100) holds no matter whether the mass hierarchy is 
normal or inverted. This is a very distinctive feature of the RG running of the leptonic mixing 
angles in the seesaw model, since such enhancement factors do not exist in the quark sector. 
Due to the threshold effects, large corrections to the leptonic mixing angles can be achieved 
even at a very low energy scale, which provides us with a possibility of testing the threshold 
effects experimentally. On the other hand, for flavor symmetric theories constructed at low 
energy scales, e.g., the TeV scale, the seesaw threshold effects should be properly taken into 
account although the RG running happens in a relatively short distance. 

Compared to the threshold corrections to the leptonic mixing angles, there is however no 
enhancement factor in the RG evolution of the neutrino masses. In the case of CP violation, 
the running CP-violation phases are also affected by the threshold effects. However, since 
there is no experimental information on leptonic CP violation, we will not present further 
analysis of the running CP-violating phases. In what follows, we will illustrate the threshold 
effects on the leptonic mixing angles and concentrate on the possibility of attaining flavor 
symmetric mixing patterns at a low energy scale. 



IV. NUMERICAL ANALYSIS OF THRESHOLD EFFECTS IN RGES 



In our numerical analysis, we adopt the values of the neutrino mass-squared differences 
Amfj and the leptonic mixing angles 8ij in the standard parametrization 34j from a global- 
fit of current experimental data in Ref. 35[ that we assume to be given at the energy scale 
\i = Mz, where Mz denotes the mass of the Z boson. The values of the quark and charged- 
lepton masses as well as the gauge couplings are taken from Ref. 36|. Furthermore, we take a 
representative value of the Higgs mass to be mu = 140 GeV. For an illustrative purpose, the 
right-handed neutrino masses are chosen as (Mi, M 2 , M 3 ) = (200 GeV, 500 GeV, 1000 GeV). 
It is worth mentioning that our choice of the right-handed neutrino masses corresponds to 
e ~ —6 x 10~ 3 between the scales M3 and M2 and e ~ —8 x 10 -3 between the scales M2 
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FIG. 2: The RG evolution of the leptonic mixing angles 623, 612, and #13 with solid, dashed, and 
dotted curves, respectively. The threshold of the right-handed neutrinos are marked by dashed- 
dotted vertical lines. At the scale M3 = 1 TeV, the bi-maximal mixing pattern, i.e., 612 = 823 = 7r/4 
together with #13 = 0, is assumed, while the corresponding rotation angles in O are given in 
the plot. In addition, the light neutrino mass spectrum is taken to be normal hierarchy with 
mi(M3) = 0.05 eV. Because of the lacking of experimental evidence of leptonic CP violation, we 
assume all CP-violating phases to be zero in our calculation. 



and Mi, cf. Eq. (fTTj) . 

In Fig. [21 we show one example of the threshold effects on the leptonic mixing angles. 
The rotation angles in O at the scale // = M3 have been labeled in the plot. In addition, we 
use the normal mass hierarchy, i.e., mi < rri2 < rri3 together with mi{Mz) = 0.05 eV, which 
translates into mi ~ 0.047 eV at the Mz scale. One observes that, with the threshold effects 
included, the exact bi-maximal mixing pattern (i.e., 812 = 623 = 45° and 813 = 0) at a very 
low energy scale could be compatible with experimental measurements. In particular, 6* 12 
decreases by about ten degrees from M3 to M\. This can be compared to a similar study 



of obtaining the bi-maximal mixing pattern at the GUT scale 37J. Such a significant effect 
can be understood from our analytical results. For instance, in the RG running between 
M 3 and M 2 , we have Sis 2 c 2 = sin(57°) sin(25°) cos(25°) ~ 0.3, e ~ -6 • 10~ 3 as well as 
v /m 1 m 2 / (m 2 — m{) ~ 43 (at fi = M 3 ), which leads to a correction 9' 12 — —4.7° according 
to Eq. (J23J) [or Eq. (JB9]) ]. Continuing the RG running between M 2 and Mi [see Eq. flB24|) 
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we obtain a total correction 9' 12 — —9.2° between M3 and Mi, which agrees approximately 
with the numerical value that is given by 9' 12 ~ —12.5°. In general, varying the value of the 
smallest light neutrino mass as well as changing the values of the Casas-Ibarra angles, we 
observe that the analytical results, which are only accurate for small corrections, are in very 
good agreement with the numerics. 

Furthermore, a non-trivial #13 can be expected at the Mz scale because of the threshold 
corrections, although we have assumed #13 = at high energy scales. Apart from the bi- 
maximal mixing pattern, other attractive flavor symmetric structures, e.g., the tri-bimaximal 
mixing pattern (i.e., 12 ~ 35.3°, 9 23 = 45° and 13 = 0) Q-Q can also be easily obtained 



through a proper choice of the rotation angles $j (see also Refs. 4lN49| for discussions on 
possible RG corrections to the tri-bimaximal mixing pattern in other frameworks). 

Note that, although all the three mixing angles do not stay at their initial input values, 
the change of 9 12 is generally larger than those of the other two mixing angles. In view of 
Eqs. ( !22l) -(!24l). #13 and $23 could only be amplified by the inverse of Am^ = ml — m\ or 
Amjj = m§ — mf, but not the inverse of Am^ = m 2 — mf, since these angles are located 
in the third column of the standard parametrization. As we have shown in the previous 
section, the boosting effect of Am^ is much more significant than those of Am^ and Am^. 
Therefore, 612 receives a more significant RG correction. In principle, if the mass spectrum 
of the light neutrinos is highly degenerate, visible changes of all the leptonic mixing angles 
can be expected. 

It should be mentioned that, in the supersymmetric framework, the RG running of lep- 
tonic mixing angles may also be enhanced by a large value of tan/3. However, it has been 



shown that, in t 



re MSSM, 9\2 generally tends to increase when running down from high 



50 



51], which means that many flavor symmetric mixing patterns are unfa- 



energy scales 
vorable. 

For completeness, we also show the RG running of the masses of the light neutrinos in 
Fig. [3j One reads directly from the figure that there is no significant threshold corrections 
to the neutrino masses, which is consistent with our analytical results. Note that we do 
not include CP-violating effects in our analysis. Now that they may affect the RG running 
behavior of mixing angles, a thorough study of the RG running CP-violating phases is 
worthy, but out of the scope of this work, and hence will be elaborated elsewhere. 

As a result of the freedom in O, various RG running behaviors of the leptonic mixing 
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FIG. 3: The RG evolution of the masses of the light neutrinos mi, 7712, and 771,3 with solid, dashed, 
and dotted curves, respectively. The corresponding input parameters are the same as those in 



angles become viable for different choices of 0j. Once a certain mixing pattern is fixed at 
a high energy scale, i.e. \i = M 3 , one can tune the rotation angles 0j so as to achieve the 
experimentally measured data at low energy scales. Concretely, we show an example in 
Fig. IH in which the bi- maximal mixing pattern of U is fixed at the scale \x = M 3 , while 
the parameter spaces of 0j at la, 2cr, and 3cr CL. are shown in the plot. The colored areas 
correspond to the choice mi(M 3 ) = 0.05 eV, while the curves mi(M 3 ) = 0.03 eV. As we 
have estimated before, in the case of mi = 0.05 eV, significant corrections could be received 
for 612, and therefore, the bi- maximal mixing pattern can be easily fetched at higher energy 
scales. For the case of a smaller mi, the parameter spaces of 0j shrink a lot, and a larger 
value of 0i, e.g., 0i ~ 90°, is more favorable. In particular, 0i ~ is ruled out, which can 
be seen from Eq. ( 123T) where the threshold correction to #12 is suppressed for small 0i. 

We remark that, although we only consider the normal hierarchy of light neutrino masses 
in our numerical analysis, the main results are unchanged if the inverted mass hierarchy 
is assumed, i.e., m 3 < mj < m 2 . The reason is simply that the RG running of #12 is 
mainly affected by the inverse of Am^, which does not flip its sign for different neutrino 
mass hierarchies. In the limit of a highly degenerate spectrum of light neutrino masses, 
i.e., mi ~ m 2 — m 3 , sizable enhancements of RG running can also be acquired for both 



Fig. El 
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FIG. 4: Parameter spaces of #j relating the bi-maximal mixing pattern at the scale fi = M3 to 
the experimentally measured leptonic mixing angles at la (red area, solid curve), 2a (orange area, 
dashed curve), and 3a (yellow area, dotted curve) C.L. The colored areas correspond to the choice 
m\{M^) = 0.05 eV, while the curves m\(M$) = 0.03 eV. The black diamonds represent the input 
values of in Figs. [2] and El 

#13 and 623, of which the light neutrino mass hierarchies may play a key role in the RG 
running. Furthermore, in our numerical calculations, we have fixed the mass spectrum of 
the right-handed neutrinos, in which the interval between two such neutrinos is around a 
few hundred GeV. In general, a stronger hierarchy among the masses of the right-handed 
neutrinos can lead to more sizable RG corrections to the neutrino parameters, whereas a 
weakly hierarchical mass spectrum of the right-handed neutrinos would suppress the seesaw 
threshold effects. 

V. SUMMARY AND CONCLUSIONS 

In this work, we have studied the RG running of neutrino parameters in the low-scale 
type-I seesaw model. We have shown that, in the RG evolution of the leptonic mixing angles 
between the seesaw thresholds, significant radiative corrections could be obtained even for a 
short distance of RG running. In particular, we have presented analytical formulas for the 
RG corrections to the neutrino parameters in crossing the seesaw thresholds, which clearly 
indicate that sizable threshold corrections to the leptonic mixing matrix may occur due to 
the mismatch between different contributions to the mass matrix of the light neutrinos. In 
addition, we have demonstrated that, if the light neutrino mass spectrum is nearly degener- 
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ate, there exist large enhancement factors boosting the RG running of the leptonic mixing 
angles remarkably. A numerical example has also been given to show that, in the presence 
of low-scale right-handed neutrinos, the bi-maximal mixing pattern at TeV scale is fully 
compatible with the current measured leptonic mixing angles. This peculiar feature turns 
out to be very useful for the construction of a low energy scale flavor symmetric theory, and 
in particular, it provides us with the possibility of testing such a flavor symmetry at running 
and forthcoming experiments (such as the LHC). For the sake of keeping the seesaw scale 
within the reach of these experiments, the neutrino Yukawa coupling tends to be chosen 
very small in order to stabilize the mass scale of the light neutrinos. One might think that 
this prevents the low-scale type-I seesaw model from being testable, since the predominant 
interaction between the right-handed neutrinos and the SM fields is through the Yukawa 
coupling. However, this is not necessarily true in the presence of additional interactions. 
For example, if the right-handed neutrinos and other SM fields are charged under an addi- 
tional gauge group, broken not far above the TeV scale, such a new gauge interaction could 
give rise to observable effects. Note that, although we mainly concentrate on the type-I 
seesaw model throughout this work, our main conclusions can be easily extended to various 
low-scale seesaw models. For example, in the type-I+II seesaw framework with three right- 
handed neutrinos and a heavy Higgs triplet field being introduced, both the right-handed 
neutrinos and the Higgs triplet may contribute to the masses of the light neutrinos, and 
the mismatch between different contributions would result in similar threshold corrections 
to the leptonic mixing angles. 
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Appendix A: Derivation of the analytical formulas for the corrected leptonic mixing 
matrix and the light neutrino masses 



In this appendix, we derive the first-order threshold corrections to the leptonic mixing 
matrix and the light neutrino masses. In general, the light neutrino mass matrix is complex 
and symmetric, i.e., m v = m%. In order to derive the series expansions, we write the 
perturbed mass matrix as 

m v = m u + ml, (Al) 

where m' v corresponds to a correction to m v . The corrected leptonic mixing matrix is given 
by U = U + U' with U' being the correction of U. In analogy, one can define the real and 
diagonal matrix D = D + D', i.e., D = diag(m 1; m 2 > ^3), and similarly for the matrices D 
and D'. 

By definition, we have 

UDU T = fh v , UDU T = m v , (A2) 

and to lowest order in the corrections U' and D', we find that 

m ' v = mu - mu = UDU T - UDU T ~ UDU' T + U'DU T + UD'U T . (A3) 

Multiplying with U> from the left and U* from the right to the above equation, we obtain 

Xflm'JU* = DU' T U* + tfU'D + D'. (A4) 

Now, we define an anti-Hermitian matrix T and a symmetric matrix A as 

T = rfU', (A5) 
A = U ] m' v U\ (A6) 

Here one can easily prove that = — T, and the diagonal elements of T are purely imaginary, 
i.e., Ta = —T* t . Using Eqs. (1A4j) - (lA6[) . A can be computed as 

A = —DT* +TD + D'. (A7) 

In order to obtain the corrections to the mass eigenvalues, we investigate the diagonal part 
of this matrix equation, which indicates that (no summation over repeated indices) 

An = -rriiT^ + T^m; + m-, (A8) 
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and the real part gives the corrections to the neutrino masses, i.e., 

m'i = Re(A ii ). (A9) 

In addition, the imaginary part is 

*)^. (A!0) 

Next, to obtain the corrections to the leptonic mixing matrix, we study the non-diagonal 
part of Eq. flA7j) . which is 

Aij = -rriiT*! + T^rrij, i ^ j. (All) 
Now, the real and imaginary parts of Eq. flAllj) yield 

= (A13) 



respectively. Therefore, all the elements of T can be obtained by using Eqs. OAlOp . (IA12[) 
and ( 1A13I) . and U' can be calculated according to Eq. ( 1A5I) . 

Note that, if m' u = £,m u , where ^ is a constant, then D = (1 + £)D and T = 0, and 
thus, U — U. This means that, if m' v is proportional to m v , U will not be affected. As a 
special case, consider CP conservation, when the original mixing matrix U as well as the 
perturbation m' v are real. Then, A* = A, and therefore, one has 

m'i = An, (A14) 
Tn = 0, (A15) 

Tij = — — — , (A16) 
rrij — mi 

As one application, we consider the threshold corrections in the main context. Below the 
scale fi = M n , one should take [cf. Eq. (TTJJ] 

m' v = (b-l)m u + ev 2 K {n) , (A17) 

as the perturbation to m v . Using the matching condition at [i = M n , can be expressed 
in terms of Y v as 

v (n) v (n)T 

«W = , (A18) 
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where y„ is the n-th column of Y v . Here we will make use of the Casas-Ibarra parametriza- 
tion of Y v [cf., Eq. (ITS]) ], which allows us to express the corrections to U in terms of the 
rotation angles Hi and the low-energy parameters. For a general U and Yj,, one can calculate 
the matrix A = A^ n \ giving the contribution of for u below the threshold M n , 2 



\\f = {U ] m' y U*).. = (6 - l)<Jij77ii + eO in O jny /m~W j 



(A19) 



Inserting the above equation into Eq. flA9j) . one directly arrives at the neutrino masses at 
the scale u = M n _i, i.e., 



m- 



[b + eRe(OD] 



rrii. 



(A20) 



Also, the components of the matrix T can be calculated as 



in) 



-lm(0 
2 v 



2 ) 

inn 



rpyi 



(n) 



Re(O in O jn ) .Jm(O in O jn ) 



rrij — rrii 



rrii + i 



(A21) 
(A22) 



Note that Oij may, in principle, have very large absolute values if they are complex. If O is 
real, or purely imaginary, 3 the corrected mixing matrix at the lower scale \x = M n _\ reads 



Uai = U aj + [UTW] aj = U aj + e 



y/rrijmk 
rrij — rrik 



(A23) 



Appendix B: Corrections to the leptonic mixing angles and the light neutrino 



masses 



In the following, we use the analytical results obtained in Appendix |X] to derive explicit 
expressions for the threshold corrections to the leptonic mixing angles and the light neutrino 
masses. Here, for simplicity, we assume CP is preserved, i.e., all the CP-violating phases are 
taken to be zero. Since the perturbation is linear in k, one should, for \i between the n-th 
and (n — l)-th thresholds, simply add the contributions from all the thresholds above u. 



2 If all the right-handed neutrinos are located at the same energy scale, one would effectively obtain a 
summation over the index n, giving the result Ay = (b — 1 + £)5ij7rii = aSijrrii, since OO t = 1. 

3 That is with "low-energy" but without "high-energy" CP-violation. 
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1. Threshold corrections between M 3 and M 2 



a. Correction to 9i 3 

Using Eq. 0A23|) with n = 3, the corrected quantity to S13 is given by 



(3) \ Jm 3 m] t 
s i3 = s i3 - S13 = U e3 ~U e3 = } U ek T^' = e y — o k3 o 33 u ek 

k^3 kj£3 d K 

Jm x m 3 X3 33 Jm 2 m 3 23 33 

= e- C13C12 + e- C13S12 

m 3 — mi m 3 — mi 



. s 2 ci2^mi Sic 2 si2y/m 2 . 

£cic 2 ci 3 < s /m 3 1 . (Bl) 

m 3 — mi m 3 — m 2 



Taylor expanding this equation around s 13 = s 13 , one obtains 



•5'', 



6 13 = arcsin( Sl3 + s' 13 ) = 6 13 + + O [(4) 2 ] 

Cl3 



~ 013 + ECiC 2 ^m 3 1 . (B2) 

y m3 — ?7i2 m 3 — mi J 
In the nearly degenerate case, i.e., mi ~ m2 ~ 7713, the above formula approximates to 



or a a C!C 2 (S2C12 + S1C2S12) v/mlmi 

"l3 = "l3 — #13 — ^ • (odj 

m3 — mi 

One can observe that the correction to 8i 3 is actually independent of S13. Furthermore, note 
that this expression is proportional to y/m\m 3 / (777.3 — mi). Thus, the correction cannot be 
very large, since there is no l/(m 2 — mi) enhancement. 

b. Correction to 9i 2 

Similarly, using Eq. (1A23|) . the threshold correction to U e2 can be obtained as 



/ ~ \— ^ (3) «, yjm 2 m k 

u e2 = u e2 — u e2 = y u e kT k2 = e y Ok 3 o 23 u e k 

z — ' ' m 2 — m k 

— /s 2 Ci2Ci 3v /mT Cic 2 si 3v /m3"\ 

= esiC 2 Jm 2 - - . (B4) 

\ m 2 — mi m 3 — m 2 J 

To obtain the correction to si 2 , we note that, in the standard parametrization, 

S U = — • (B5) 

C13 



19 



Expanding in the small correction s' 13 , we have 



+ (B6) 



Cl3 C13 V CJ 3 

and therefore, we find that 

Sl2 = {Ue2 + U e2 ) H g - S 12 + S12 5 1 ' ( B7 ) 

Cl3 V C 13 / C l3 C13 

where the corrections s',, and t/'n have been calculated above. For ~ 0, one obtains 



Jvn\vn 2 

S12 = S 12 - S12 ^ £SiS2C 2 Ci2 . (B8 

?7l2 — ?7li 



In addition, Taylor expanding around s 12 = s i2 gives 

0' 12 = 12 - 12 ~ es lS2 c 2 ^^ , (B9) 

m 2 — mi 

where the right-hand side is independent of 8i 2 . 
c. Correction to 623 

Finally, using Eq. (IA23j) . the threshold correction to is given by 



= Up - Up = £ U, k TjS = e £ ^^O k3 33 U, k 



m 3 — m k 

fc^3 fc^3 J 



, — f s 2 U IMl ^/rrH s l c 2 U il2 ^/m 2 ~\ m . 

= ecxOiJmi — 1 — . (BIO) 

\ m 3 — rti\ m 3 — m 2 J 

As for the corrections to S23, we note that, in the standard parametrization, 

§23 = 5* (BH) 
Cl3 

and thus, one obtains 

§23 = (U, 3 + C/; 3 ) — ( 1 + ^1 * *23 + ^23^ + ^, (B12) 
c 13 V c 13 / C 13 C 13 

where s' 13 and U' 3 have been calculated above. For s i3 ~ 0, one obtains 



'23 - s 23 - «23 - tc\C2C 2 -iyfmz s 2 s 12 h s x c 2 c X2 , (B13) 

m 3 — rrii m 3 — "Hi 2 



giving that 



#23 = #23 - 23 ~ ecxc 2 ^fm3 -s 2 syi h s x c 2 c 12 . (B14) 

\ m 3 — mi m 3 — m 2 J 
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One finds that 9' 23 is independent of #23- In the nearly degenerate limit, the above formula 
approximates to 

Jm\m 3 

6> 23 = £CiC 2 (S1C2C12 - S 2 S 12 ) • (B15) 

m 3 — mi 

Similar to the corrections to #13, there is no strong enhancement factor. 



d. Corrections to the light neutrino masses 

According to Eq. (IA20j) . the explicit expressions for the threshold corrections to the light 
neutrino masses are 



rn'i 


= mi 


— mi = 


[(»- 


i)- 


1- esl] mi, 


m' 2 


= m 2 


- m 2 = 


[(»- 


i)- 


1- es\d£[ m 2 


m' 3 


= rh 3 


- m 3 = 


[(»- 


i)- 


1- ec\c^\ m 3 



(B16) 
(B17) 
(B18) 



2. Threshold corrections between M 2 and M\ 



One may repeat the above analysis with n — 2, and calculate the threshold corrections 
to the neutrino parameters between M 2 and M\. Explicitly, for S13, one obtains 

/ n 1 ( C2S3C12V 7 ™! . 22 S 12 y/mi\ , . 

s 13 = e0 32 c 13 ^/mi 1 , (B19) 

\ m 3 — mi m3 — m 2 / 

and henceforth, 

l3 = e0 32 ^m7 3 - + , (B20) 

y m 3 — mi m 3 — m 2 / 

where O22 = QlC 3 — S\S 2 s 3 and O32 = — S1C3 — C1S2S3. In the nearly degenerate limit, the 

above expression approximates to 



ff ~ e0 32 v/mim3 ( C2S3 ci2 + 22 s 12 ) . (B21) 
m 3 — mi 



Then, for #12, one obtains 



„, n , — ( c 2 s 3 c 12 c 13 ^/m^ 32 s 13 ^/Wi\ 

u e 2 = t0 22 Jm 2 . (B22) 

\ m 2 — mi m 3 — m 2 



Ignoring the small parameter s 13 , we arrive at 



/ Jmim 2 

s 12 ~ ec 2 s 3 (cic 3 - S1S2S3) C12 , (B23) 

m 2 — mi 
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which corresponds to 



9 12 ~ ec 2 s 3 (cic 3 - S1S2S3) • (B24) 

m 2 — mi 



Next, for # 23 , the correction to reads 



p V m 3 — mi m3 — m 2 



which, for si 3 ~ and in the nearly degenerate case, gives 



and 



s' 23 ~ e0 32 (-C2S3S12 + 22 c 12 ) c 23 ^ mim3 (B26) 

m 3 — mi 



£ 23 ~ e0 32 (-c 2S3 si 2 + 22 ci 2 ) Vmim3 . (B27) 

m3 — mi 



Finally, the threshold corrections to the light neutrino masses are given by 



m [ = [(b - 1) + ec 2 2 sl] mi, (B28) 
J 2 = [{b - 1) + e (cic 3 - s lS2 s 3 ) 2 ] m 2 , (B29) 



rn 



m' 3 = [(b- 1) + e (sic 3 + cis 2 s 3 ) 2 ] m 3 . (B30) 
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